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Abstract. We have solved the one-dimensional stationary two-fluid hydrodynamic equations for post-shock flows 
on accreting magnetic white dwarfs simultaneous with the fully frequency and angle-dependent radiative transfer 
for cyclotron radiation and bremsstrahlung. Magnetic field strengths B = 10 to 100 MG are considered. At given 
B, this theory relates the properties of the emission region to a single physical parameter, the mass flow density 
(or accretion rate per unit area) rh. We present the normalized temperature profiles and fit formulae for the peak 
electron temperature, the geometrical shock height, and the column density of the post-shock flow. The results 
apply to pillbox-shaped emission regions. With a first-order temperature correction they can also be used for 
narrower columns provided they are not too tall. 
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1. Introduction 

The thermal structure of the accretion columns on accret- 
ing magnetic white dwarfs can be derived analytically for 
a single-particle fluid and sufficiently simple assumptions 
on the radiative cooling (Aizu 1973, Chevalier & Imamura 
1982, Wu et al. 1994). For the more general case of the 
optically thick frequency and angle-dependent radiative 
transfer in a two-fluid plasma, the coupled hydrodynamic 
and radiative transfer equations have to be solved numer- 
ically (Woelk & Beuermann 1996, henceforth WB96). In 
this paper, we present results which are improved and ex- 
panded over those of WB96. We obtain the temperature 
and density profiles for plane-parallel post-shock cooling- 
flows and derive fit formulae for the peak electron tem- 
perature T max , the column density x s , and the geometrical 
shock height hgh as functions of the magnetic field strength 
B and the mass flow density (accretion rate per unit area) 
rh. For low rh and high £?, we show that the shock solution 
merges into the non-hydrodynamic bombardment solution 
for an atmosphere which is heated by a stream of fast ions 
and cools by cyclotron radiation (Woelk & Beuermann 
1992, 1993, henceforth WB92, WB93). 

Our treatment of radiation-hydrodynamics is one- 
dimensional and stationary. The one-dimensionality im- 
plies that our solutions are strictly applicable only to 
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pillbox-shaped emision regions with a width D 3> /i s h 
and a stand-off distance h s h "C R w d, where i? w( j is the 
white dwarf radius. The stationarity implies that our solu- 
tions describe the mean properties of the shocks and that 
aspects like rapid fluctuations in the mass flow density 
and the stability against shock oscillations are left aside. 
Shock oscillations have been treated by a number of au- 
thors (Imamura et al. 1996, Saxton & Wu 1999, and refer- 
ences therein) and generally suggest that cyclotron cooling 
stabilizes the flow and bremsstrahlung cooling destabilizes 
it. Observationally, optical oscillations have been found in 
a few polars, while the search for hard X-ray oscillations 
has so far yielded only upper limits (Larsson 1992, Wolff 
et al. 1999, Imamura et al. 2000, and references therein). 

2. Two-fluid radiation hydrodynamics 

2.1. General approach 

We solve the stationary, one-dimensional, two-fluid hydro- 
dynamic equations simultaneous with the frequency and 
angle-dependent radiative transfer, closely following the 
approach of WB96. We deviate from WB96 in the treat- 
ment of the shock itself. Instead of integrating the flow 
through the shock with an artificial viscosity, we adopt 
the presence of a strong ion shock and start the integration 
with adopted values of the ion and electron temperatures 
(see below). Of course, the solution now fails to reproduce 
the rapid rise in ion temperature across the shock, but 
otherwise the results are practically identical except for 
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small differences at low velocities where the flow connects 
to the atmosphere of the star and large gradients revive 
the viscous terms again. The set of differential equations 
then reads (compare Eqs. (1) to (4), (7), and (8) of WB96) 



for emission regions of finite lateral extent D (Fig. 1) is 
discussed in Sects. 3.2 and 3.2 below. 



pv = —m 

p v ^r + -7-( p i + p c) = -9 

dx dx 

d£ c , D . e , dp dF rad 
pv -. v (P c + c c )— = A ci - p — — 

da: dec dx 



(1) 
(2) 

(3) 

(4) 



pcostf^^ = + a] I,A$) - K v {$)B v (T e ) ~ a J, (5) 



dx 

oo +1 

F rad = 2tt J J J v (i?)cos0d(cosi?)di/ 

-1 



(6) 



where p is the mass density, v the velocity, m the mass 
flow density, Pi and P c the ion and electron pressures, E\ 
and £ c ion and electron internal energy densities, I v {"&) the 
specific intensity of the radiation field at frequency v and 
angle J v the angle-averaged intensity, the angle- 

dependent absorption coefficient, a the Thomson scatter- 
ing coefficient, B V (T C ) the Planck function at electron tem- 
perature T c , and F ra d the total radial energy flux in the 
radiation field. We choose the downstream column den- 
sity x as the independent variable rather than the radial 
coordinate or the geometrical height h. Following WB96, 
we neglect the effects of radiation pressure and thermal 
conduction. 

The connecting link between the hydrodynamics (Eqs. 
1 to 4) and the radiative transfer (Eqs. 5 and 6) is 
Fi-ad: the electron gas cools by radiation and is heated by 
Coulomb interactions with the ions, described by the non- 
relativistic electron ion energy exchange rate A e ; (Spitzer 
1956, see also WB96, their Eqs. 5 and 6). The fully angle 
and frequency-dependent radiative transfer accounts for 
cyclotron absorption, free-free absorption, and coherent 
electron scattering. Our emphasis is on the largely correct 
treatment of the cyclotron spectra and we accept inaccu- 
racies of the hard X-ray spectra caused by the neglect of 
Compton scattering. This still rather general treatment 
ensures that our results are relevant for a wide range of m 
including the low-m regime where radiative losses by op- 
tically thick cyclotron radiation dominate. The cyclotron 
absorption coefficients used here are the added coefficients 
for the ordinary and the extra ordi nary rays (WB92). This 
limitation is dropped in Sect. 3.1, below. 

We use a Rybicki code for the LTE radiative trans- 
fer and integrate the set of equations implicitly, using a 
Newton scheme to iterate between hydrodynamics and ra- 
diation transport. For more details see WB96. Our solu- 
tion is strictly valid only for an infinite plane parallel layer. 
A first-order correction to the peak electron temperature 



Eq. (2) accounts for post-shock acceleration and heat- 
ing of the flow by the constant gravity term g = 
GM wd /R wd 2 . Within our one-dimensional approximation 
which disregards the convergence of the polar field lines, 
considering the variation of gravity with radius would 
not be appropriate. Our approach is, therefore, limited 
to stand-off distances of the shock /i s h < Ci? w d- Settling so- 
lutions with h s h <; i? w d are not considered. 

As in WB96, we assume that the pre-shock flow is fully 
ionized, but cold. Soft X-rays will photoionize the infalling 
matter and create a Stromgren region with a temperature 
typical of planetary nebulae, but for our purposes this is 
cold. Heating of the pre-shock electrons by thermal con- 
duction may be more important. Equilibrium between dif- 
fusion and convection defines an electron precursor with 
a radial extent A pre — 4 x lO~ 15 Tc/ s 2 jm cm, where T c s 
is the electron temperature at the shock in K (Imamura 
et al. 1987) and m is in gcm~ 2 s _1 . Near the one-fluid 
limit, electron and ion shock temperatures are similar, 
T CiS ~ T; iB) and the precursor extends to A prc ~ 0.09 /i s h- 
In a cyclotron-dominated plane-parallel flow, however, two 
effects cause the precursor to be less important: (i) the 
electrons never reach the peak temperature expected from 
one-fluid theory and (ii) the optically thick radiative trans- 
fer in the plane-parallel geometry sets up a radial temper- 
ature gradient which further depresses the electron tem- 
perature at the shock. In this paper, we do not consider 
thermal conduction, neglect the presence of the electron 
precursor, and opt to set T @tS = 0. 

At x = 0, we adopt the Rankinc-Hugoniot jump con- 
ditions for a gas with adiabatic index 5/3, i.e. we set the 
post-shock density to 4p , the bulk velocity to v a /4, and 
the pressure to (3/4)p v Q 2 , with p Q and v Q the density 
and bulk velocity in the pre-shock flow. With T e<s = 0, the 
ion shock temperature is 



T iB = 3^im u /(16k)Uo. 



(7) 



We use u = (2GA/ wd /i?wd) 1/2 with Nauenberg's (1972) 
relation between mass and radius of the white dwarf. p\ is 
the molecular weight of the ions, m u the mass unit, and k 
the Boltzmann constant. 

All numerical calculations are performed for a hydro- 
gen plasma with p = 0.5 and p\ = fi e = pc% = 1, where p e 
is the number of nucleons per electron, \i is the molecular 
weight of all particles, and pz the molecular weight of the 
ions weighted with Z\. We include the molecular weight 
dependence in our equations in order to allow conversion 
to other compositions, e.g., a fully ionized plasma of solar 
composition with p = 0.617, p\ = 1.297, p c — 1.176, and 
pz = 0.927. 
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Fig. 1. Schematics of the emission region. The region is 
bounded at the top by the shock front and at the bottom 
by the white dwarf. 

2.2. Bremsstrahlung and cyclotron emissivities 

The frequency-integrated volume emissivity for brems- 
strahlung is 



£brc 



c 2 T 



V2, 



C2 



1/2 



^ nkZk ^(km3)i/2 Mz 



pvy/2( 8 ) 



where C2 = 1.43 1CT 27 cgs-units, n c is the 
number density of the electrons, the 
number density of the ions of charge Z^e, 
P is the gas pressure, and p the mass density. The 
Thomson scattering optical depth of a bremsstrahlung 
dominated flow parallel to the flow is of order unity, 
implying that bremsstrahlung is essentially optically thin 
(Aizu 1973). 

We use the cyclotron absorption coefficients for the 
ordinary and extraordinary rays given by Chanmugam & 
Dulk (1981), Thompson & Cawthorne (1987), and WB92. 
The total cyclotron emissivity of non-relativistic electrons 
integrated over wavelength and solid angle is 



e cyc oc n Q T B 2 oc (fi/fi e ) B 2 P. 



(9) 



In the columns considered, cyclotron radiation is optically 
thin in the higher harmonics, but is always optically thick 
in the first few harmonics, and the temperature distribu- 
tion in a cyclotron-dominated emission region can prop- 
erly be calculated only by solving the coupled radiation- 
hydrodynamic equations. 

2.3. Geometry of the emission region 

Fig. 1 shows the schematic of an emission region with finite 
lateral extent D. The shock is located at h = ft. s h above 
the white dwarf surface. The downstream column density 
is x = at the shock and and x — x s at the surface of the 
star, with x and h being related by da; = —pdh. The grav- 
ity vector g and the magnetic field vector B are taken par- 
allel to the flow lines. The radiation- hydrodynamic equa- 
tions are solved for layers of infinite D to yield the run 
of electron temperature and mass density, T e (x) and p(x). 
These profiles are later employed to calculate the outgoing 
spectra for emission regions with finite D by ray tracing, 



i.e. by adding the contributions from an appropriate num- 
ber of rays (Fig. 1 and Sect. 3.1). This procedure is not 



self-consistent if optically thick radiative losses occur from 
the sides of the column. An appropriate first-order correc- 
tion to the temperature structure derived for the infinite 



layer is discussed in Sects. 3.2 and 3.3, below. The treat 



mcnt of really tall columns requires a different approach 
which specifically allows for the emission from the sides of 
the column (Wu ct al. 1994). 

Radiation intercepted by the white dwarf is either re- 
flected or absorbed and reemitted by its locally heated 
atmosphere. We assume coherent scattering of hard X- 
rays using the frequency-dependent reflection albedo A v 
of van Teeseling et al. (1994). The fraction 1 — A v of the 
energy is re-emitted in the UV and soft X-ray regime and 
is not considered in this paper. 

2.4. Limiting cases 

Here, we consider simple limiting cases which can, in part, 
be solved analytically. Below, we shall discuss our numer- 
ical results in terms of these limiting solutions. The high 
to, low B limit is the bremsstrahlung-dominated one-fluid 
solution. In the opposite limit of low to, high B one enters 
the non- hydrodynamic regime (Lamb & Masters 1977). 
Here, the bombardment solution of a static atmosphere 
heated by a stream of fast ions and cooling by cyclotron 
emission is an appropriate approximation (Kuijpers &i 
Pringle 1982, WB92, WB93). 

The one-dimensional, one-fluid hydrodynamic equa- 
tions with simple terms for optically thin cooling can 
be solved analytically (Aizu 1973, Chevalier & Imamura 
1982). Integration of Eq. (2) with Eq. (1), P = P + P e , 
and g — yields P — m(v — v) which allows to ex- 
press the emissivities of Sect. 2.2 as ebrcms oc to 2 f(u) and 
e cyc oc to g(u), with f and g being functions of the flow ve- 
locity v and with additional dependencies on the /i's and 
B contained in the proportionality factors. Integration of 
the energy equation over v yields expressions for the col- 
umn density a; s and the geometrical shock height h s h which 
reflect the parameter dependence of e, 
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constant 
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ecyc- 



In the flows considered here, bremsstrahlung is close to 
optically thin and the analytical solution is quantitatively 
corroborated by our numerical results. Cyclotron emis- 
sion, on the other hand, is optically thick in the lower har- 
monics which reduces the effective emissivity and inflates 
the emission region. While the TO-dependence of cyclotron- 
dominated columns in Eqs. ( |To| ) and (11) is recovered in 



our numerical calculations, not surprisingly, the numer- 
ically derived sizes of x s and h s h are much larger than 
those predicted by the (unquoted) proportionality factors 
for cyclotron cooling in Eqs. ( |l0| ) and ([Tl|). 
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Fig. 2. Temperature profiles for the ions (dashed curves) and electrons (solid curves) as functions of column density 
x for M wd = 0.6 M Q and field strengths of 10 MG (left panel), 30 MG (center panel) and 100 MG (right panel). The 
individual curves are for mass flow densities rh = 100, 10, 1, 10 _1 and 10 _2 gcm _2 s _1 (from top). In the left panel, 
the curves for 100 and 10 gcm _2 s _1 are indistinguishable. In the right panel, the bottom curve is for 10 _1 gcm" 2 s" . 
From Eqs. (|l|) and @, z s ,brcms = 0.783gcm~ 2 (log 2 s ,brems = -0.106) and T maXjbrems = T ijS /2. 



2.5. Bremsstrahlung-dominated shock solution 



2.6. Cyclotron-dominated bombardment solution 



For a strong shock in a one-fluid plasma with adiabatic 
index 5/3, the normalized post-shock velocity u) — 4u/d 
varies between 1 and 0. The column density x measured 
from the shock is related to u by (Aizu 1973, Chevalier & 
Imamura 1982) 



ci« 2 



V3- 



1 



■y/4w - 1 



cos 



1 



3 2 

The total column density and shock height are given by 

(13) 
(14) 



The bombardment solution involves by nature a two-fluid 
approach. WB92 solved this case using a Fokker-Planck 
formalism to calculate the stopping length of the ions and 
a Feautrier code for the radiative transfer. WB93 (their 
Eqs. 8, 9) provided power law fits to their numerical re- 
sults for the column density and the peak electron tem- 
.(12j))erature. Since the ions are slowed down by collisions with 
atmospheric electrons, a factor /i c appears in x s : 



^s,brcms 
^sh,brems 



-2 



(\/3-7T/3) g cm 

(39\/3 — 207r^ /(48 rh) 
1 a „„\f ,, ,,_/,, 1/2-1 



cm, 



of 



where c\ = (k 1 / 2 mu /z /4c2)(/x e /zz/A' 1 ^ 2 ) with values 
6.22 10 -18 cgs for pure hydrogen and 6.10 10 -18 cgs for 
solar composition. The temperature profile follows from 
pressure equilibrium P — pv(v — v), the equation of state 
for the ideal gas, and Eq. (1) as 



T 



T, 



max.brcms 



with 

^max,brcms 



3 /zm u / (16 k) vl = (p/ p{) T ijS 



(15) 



(16) 



Our two-fluid calculations for high rh, low B reproduce 
the temperature profile T(x) given by Eq. ([l5|) with ( |l2| ) 
and ( |l6| ) to better than 1 % of T max ,brems) except for the 
initial equilibration layer which is infinitely thin in the 
analytic calculation and has a finite thickness with rising 
electron temperature in our calculations. As an aside, we 
note that T/T maXi b r oms — (1 — x/x s ) ^ 9 with an rms error 
of less than 1%. 



z s ,bo,„b - 3.94 x 10- 2 M e (mS7 2 - 6 ) 30 M4 d 72 



max, bomb — 1-28 X 10 (mB-j 



2.6\0.42 i f 0.66 
) m wd 



K. 



;cn<lt) 
(18) 



Here, m is in gcm _2 s _1 , By is in units of 10 7 G, and M w d 
is in solar masses. These fits are very close to the quasi- 
analytical expressions of Eqs. (5) and (6) of WB93[]. 

With increasing rh, a shock develops which is initially 
cyclotron-dominated and ultimately bremsstrahlung- 
dominated. Since T roaXi bomb reaches T maX! b rcms at some 
intermediate rh, we expect a smooth transition in peak 
temperature between these cases. The situation is quite 
different for x s , however. At the rh where T max bomb equals 
Tmax.brems, ^s.bomb and :c s ,brems differ by more than two 
orders of magnitude. The run of x s (rh) between these two 
limiting cases can be determined only with a radiation- 
hydrodynamical approach. 

The bombardment solution does not predict the geo- 
metrical scale height of the heated atmosphere which we 
expect to lie between that of a corona with an external 
pressure P — and that of a layer compressed by the ram 
pressure P — p v 2 . 



Note the misprint in Eq. (6) of WB93 which should read 
B 085 instead of B 185 . 
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2. 7. Parameterization of the results 

In the bombardment solution (Eqs. 17, 18), the depen- 
dence of x s and T m ax on mB~ 26 is obtained from the equi- 
librium between the energy gain by accretion, .F acc oc m, 
and the energy loss by optically thick cyclotron radiation, 
F cyc cx T max u>l, where = m*ui c is the high-frequency 
cutoff of the cyclotron spectrum and uj c oc B the cyclotron 
frequency. We determine the limiting harmonic num- 
ber from the cyclotron calculations of Chanmugam 
& Langer (1991; their Fig. 5) as m» =4.43A|T| with 
a ~ 0.12 and b ~ 0.40. This approximation is valid near 
depth parameters A4 = A/10 4 ~ 1 and temperatures 
Ts = T , max /(10 8 K) ~ 1 and is more adequate for the 
cyclotron-dominated emission regions on polars than the 
frequently quoted formula of Wada et al. (1980). Replacing 
a;* with A cx x s /B in F cyc and equating the accretion and 
radiative energy fluxes yields the result that a power of 
Tmax is proportional to m£? 3 ( a_1 ) ~ mB~ 26 . The same 
holds for x s . 

We find that the cyclotron-dominated shocks at low m 
behave similarly to bombarded atmospheres in that their 
thermal properties, too, depend on rhB~ 26 . The individ- 
ual temperature profiles T c (x) for different rh, B with the 
same rhB~ 2e coincide only in an approximate way, but 
the dependency on rhB~ 2 e holds quite well for the two 
characteristic values of each profile, T max and x s . If we 
leave the exponent a in rh B a as a fit variable, the small- 
est scatter in T max and X$ clS functions of rh B a is, in fact, 
obtained for a = —2.6 ± 0.2. 

Our model calculations cover magnetic field strengths 
B = 10 - 100 MG and mass flow densities rh = 10~ 2 - 10 2 
gcm _2 s _1 . In what follows, we present first the tempera- 
ture and density profiles along the flow lines. From these, 
we obtain T max , x s , and /i s h as the characteristic parame- 
ters of the post-shock flow which are presented in an ap- 
propriate way as functions of rhB~ 26 . 

2.8. Electron temperature profiles for an infinite layer 

Fig. H shows the temperature profiles Ti(x) and T e (x) for 
Mwd— 0.6 Mq and several rn-B combinations on a log- 
arithmic depth scale which emphasizes the initial rise of 
the profiles. These profiles display a substantial spread in 
:r s and in T max , reflecting the influence of cyclotron cool- 
ing. At 10~ 2 gcm _2 s _1 , 100 MG, cyclotron cooling has re- 
duced T max to 6% and x s to 0.3% of the respective values 
for the pure bremsstrahlung solution. We have confidence 
in our numerical results because they accurately reproduce 
the analytic bremsstrahlung solution (see above). 

Fig. U displays the normalized profiles of the elec- 
tron temperature, T c /T max vs. x/x s , for different rh,B 
combinations, covering the range from a bremsstrahlung- 
dominated flow with lOOgcm^s" 1 , 30 MG (fat solid 
curve) to 10 -2 gcm _2 s _1 , 100 MG near the non- 
hydrodynamic limit (dotted curce) . They represent an ap- 
proximate sequence in mB -2,6 , but not surprisingly, the 
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Fig. 3. Normalized electron temperature distributions for 
M w( |= 0.6 Mq and the values of B and rh given in the 
figure. 

shapes differ somewhat for different rh and B combina- 
tions with the same value of rhB~ 26 (not shown in Fig. ||). 

Equilibration between electron and ion temperatures is 
reached at column densities of ~ 10~ 3 . . . 10 _1 gcm _2 s _1 
depending on rh and B (Fig. 2). At lOOgcm^s -1 , 10 MG, 
electrons and ions equilibrate as early as ~ 0.03 x s , while 
at 10~ 2 g cm _2 s _1 , 100MG, equilibration length and x s 
are of the same order, indicating the approach to the non- 
hydrodynamic regime. A peculiar feature of the latter pro- 
file is the extended low-temperature tail which was not ad- 
equately resolved by WB96. This tail appears when equi- 
libration occurs near the temperature at which cyclotron 
cooling becomes ineffective and the density is sufficiently 
high for bremsstrahlung to take over. It is hydrodynamic 
in origin. Apart from the tail, the temperature profile at 
10 -2 gcm -2 s -1 , 100 MG is very close to that obtained by 
the non-hydrodynamic approach of WB92, WB93. The 
low-tempcrature tail is responsible for a low-temperature 
thermal emission component with kT <^ 1 keV. 

The initial rise of the individual temperature pro- 
files is similar and is very rapid following approximately 
T c cx X 35 (Fig. |). One half of T max is reached at 0.001 x s 
in the bremsstrahlung-dominated case and at 0.006 x s near 
the non-hydrodynamic limit. Further downstream the pro- 
files differ substantially. In the bremsstrahlung-dominated 
case, the peak electron temperature is reached quickly, 
while in the cyclotron-dominated flow it occurs at the 
same x at which half of the accretion energy has been ra- 
diated away. The reason is that a temperature gradient is 
needed to drive about one half of the radiative flux across 
the shock front, while the other half enters the white dwarf 
atmosphere. In the plane-parallel geometry, the optically 
thick radiative transfer requires the electron temperature 
at the shock front to stay below the peak electron temper- 
ature : T CiS < T max < T brcms>max ~ 0.5 T ii8 . This is why we 
opted to start the integration with the initial values T CtS = 
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Fig. 4. Normalized velocity w — 4 v/v for M w( j= 0.6 Mq 
and the same values of field strength and mass flow density 
rh as in Fig. |^. 

and Tj s as given by Eq. (Q). Because of the rapid initial 
rise in T(x), our results would have been practically the 
same had we set T 0:S = 0.5T max . 

To facilitate the modeling of specific geometries, we 
provide the normalized temperature and density profiles 
for a sequence of rh,B combinations in Table 1. We 
also provide best fits to functions of 

mB- 26 . 



2.9. Velocity profiles for an infinite layer 

For calculations of the bremsstrahlung emission, we need 
the profiles of the mass density which varies as p cx v~ x . 
Fig. |I] shows the normalized velocity profiles for M wc i = 
0.6 Mq and the same rh, B combinations as in Fig. ||. In 
the limit of pure bremsstrahlung cooling, the velocity pro- 
file is indistinguishable from that given by the inversion 
of Eq. (|l^). Increased cyclotron cooling causes a similar 
depression at intermediate x as seen in the temperature 
profiles. Table 1 (bottom) provides the velocity profiles in 
numerical form for the same parameters as above. 



2.10. Maximum electron temperature T max 

In what follows, each model is represented by one 
"data point". Fig. | shows T max /Ti. s vs. mBf 2 ' 6 for 
M 



0.6 Mq and B = 10 - 100 MG. The dependence of 
26 is equally well documented for M wc i = 0.8 
and 1.0 Mq, but for clarity we do not show these data. The 
0.6 Mq results can be fitted by 



wd ' 

T max on rh B^ 



1 



T„ 



1 



a T, 



max, bomb 



1 



ai T t 



max.brcms 



1/a 



(19) 



with T max ,bomb from Eq. (|18|) and T maXihlems from 
Eq. (|l6|). The exponent a measures the smoothness of 
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Fig. 5. 

tion of rh B 7 



Maximum electron temperature T max as a func- 



for M 



0.6 M Q and B = 10 - 100 MG 



(B 7 = B/10 7 G). The solid curve is the fit from Eq. (|l9|), 
the solid straight line represents the bombardment solu- 
tion from Eq. @. The dotted curve and straight line 
represent the corresponding fit to the data for 1 Mq (data 
not shown). 

the transition between cyclotron and bremsstrahlung so- 
lutions. The fits for 0.6 and 1 Mq are included in Fig. || 
as the solid and dotted curves, respectively. The fit pa- 
rameters ao,a\, and a are listed in Table 2 for all three 
white dwarf masses. The fact that a\ falls slightly short 
of 1.0 indicates that the limiting value T max = (/i//ii)Tj jS 
for large rh is not yet reached at 100gcm _2 s _1 . Even at 
this high m, radiative energy losses remove some energy 
prior to equipartition. The maximum temperatures at low 
rh remain about 10% below the temperatures predicted 
by the bombardment solution of WB92,93 (straight lines, 
ao ~ 0.9). A difference as small as this is actually remark- 
able considering the substantially different theoretical and 
numerical approaches (radiation hydrodynamics vs. radia- 
tive transfer in a static atmosphere). 

2.11. Column density x s 

The transition of x s between the bombardment and the 
bremsstrahlung solutions (Eqs. [I?] and |l^) is more com- 
plicated than that of T max . Fig.TfTshows function 
of riiBf 2 6 for M wd = 0.6 Mq and B = 10 - 100 MG. 
Again, the data points for 0.8 and 1.0 M Q are not shown 
for clarity. We fit x s by 



1 



fj, 



b p c mB 7 



•£s,br< 



1/0 



(20) 



with Xs.brcms from Eq. (|13|). The fit parameters bo and f} 
are given in Table 2 for the three values of M w( j. Again, 
the fits are shown for 0.6 and 1.0 M Q (solid and dotted 
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Table 1. (a) Top: Temperature profiles T c /T max (x/ x s ) for a white dwarf mass of 0.6 M Q , field strengths of 10, 30, 
and 100 MG, and mass flow densities between 100 and 0.01 gcm _2 s _1 . First line: Analytical solution of Eqs. flF2] ) and 

(|l5|). Subsequent lines: Profiles for parameter combinations B, rh in MG and gcm _2 s~ 1 , with B7 in units of 10 7 G. 

All profiles start at (x/x s ,T e /T max ) — (0,0), end at (1,0), and are normalized to a peak value of unity, (b) Bottom: 

Same for the normalized velocity w — Av/v . These profiles start at (x/x s ,w) = (0, 1) and end at (1,0). 

rh B mB^ 2Ai T/T max vs. x/x s 

x/x s = 1(T 4 1(T 3 0.01 0.02 0.05 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 0.95 0.98 0.99 

00 1.000 0.999 0.995 0.989 0.973 0.945 0.886 0.821 0.751 0.674 0.589 0.494 0.382 0.245 0.156 0.085 0.054 

100 10 100 0.230 0.515 0.916 0.984 0.995 0.967 0.905 0.840 0.770 0.691 0.603 0.505 0.392 0.250 0.158 0.085 0.053 

10 10 10 0.230 0.515 0.915 0.984 0.996 0.968 0.906 0.842 0.770 0.692 0.604 0.506 0.393 0.250 0.159 0.082 0.048 

1 10 1 0.219 0.499 0.898 0.977 0.998 0.973 0.914 0.850 0.782 0.710 0.629 0.536 0.426 0.281 0.176 0.085 0.047 

0.1 10 0.1 0.191 0.434 0.822 0.927 0.999 0.980 0.916 0.849 0.779 0.706 0.622 0.522 0.400 0.239 0.127 0.047 0.022 

0.01 10 0.01 0.175 0.395 0.743 0.861 0.981 0.991 0.860 0.743 0.650 0.565 0.476 0.384 0.285 0.165 0.092 0.039 0.021 

100 30 5.75 0.228 0.512 0.914 0.984 0.994 0.966 0.903 0.836 0.766 0.686 0.598 0.499 0.388 0.247 0.153 0.080 0.046 

10 30 0.58 0.224 0.504 0.909 0.981 0.994 0.958 0.888 0.813 0.738 0.657 0.569 0.473 0.361 0.225 0.134 0.065 0.035 

1 30 0.058 0.209 0.468 0.865 0.958 0.997 0.937 0.815 0.717 0.629 0.546 0.460 0.372 0.272 0.162 0.091 0.042 0.025 

0.1 30 0.0058 0.185 0.408 0.738 0.849 0.971 0.997 0.833 0.657 0.548 0.468 0.401 0.336 0.274 0.197 0.145 0.085 0.045 

0.01 30 0.00058 0.178 0.390 0.710 0.802 0.922 0.994 0.878 0.652 0.484 0.362 0.292 0.245 0.196 0.119 0.074 0.038 0.022 

100 100 0.25 0.222 0.497 0.898 0.979 0.992 0.942 0.849 0.766 0.685 0.603 0.517 0.424 0.317 0.188 0.104 0.044 0.020 

10 100 0.025 0.206 0.468 0.839 0.940 0.999 0.927 0.756 0.639 0.541 0.450 0.367 0.289 0.201 0.103 0.043 0.011 0.006 

1 100 0.0025 0.199 0.440 0.795 0.903 0.997 0.921 0.564 0.458 0.375 0.302 0.237 0.193 0.129 0.058 0.027 0.013 0.007 

0.1 100 0.00025 0.180 0.354 0.626 0.727 0.880 0.976 0.966 0.600 0.131 0.097 0.084 0.077 0.063 0.047 0.038 0.030 0.022 

rh B mBj 2 ' 6 w = 4v/v vs. x/x s 

x/x s = 10" 4 10" 3 0.01 0.02 0.05 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 0.95 0.98 0.99 

00 1.000 0.999 0.992 0.984 0.960 0.921 0.841 0.761 0.678 0.594 0.506 0.413 0.311 0.193 0.121 0.065 0.041 

100 10 100 1.000 0.999 0.992 0.985 0.961 0.920 0.840 0.760 0.676 0.593 0.504 0.411 0.310 0.191 0.119 0.064 0.039 

10 10 10 1.000 0.999 0.991 0.983 0.958 0.915 0.834 0.755 0.670 0.586 0.496 0.405 0.304 0.187 0.115 0.060 0.035 

1 10 1 1.000 0.998 0.983 0.969 0.932 0.882 0.788 0.700 0.618 0.538 0.459 0.376 0.289 0.182 0.112 0.053 0.029 

0.1 10 0.1 1.000 0.993 0.934 0.888 0.794 0.696 0.571 0.480 0.404 0.336 0.279 0.218 0.158 0.089 0.046 0.017 0.008 

0.01 10 0.01 1.000 0.989 0.892 0.811 0.635 0.440 0.266 0.193 0.150 0.116 0.090 0.066 0.045 0.024 0.013 0.005 0.003 

100 30 5.75 1.000 0.999 0.990 0.981 0.953 0.913 0.831 0.751 0.667 0.584 0.497 0.402 0.303 0.186 0.115 0.059 0.033 

10 30 0.58 1.000 0.996 0.972 0.948 0.902 0.947 0.758 0.676 0.597 0.517 0.437 0.352 0.264 0.159 0.094 0.044 0.022 

1 30 0.058 1.000 0.990 0.914 0.852 0.726 0.614 0.499 0.425 0.363 0.307 0.254 0.199 0.144 0.084 0.047 0.019 0.009 

0.1 30 0.0058 1.000 0.990 0.902 0.829 0.664 0.480 0.286 0.206 0.167 0.140 0.117 0.097 0.078 0.055 0.040 0.021 0.011 

0.01 30 0.00058 0.999 0.982 0.895 0.824 0.653 0.435 0.148 0.087 0.062 0.045 0.036 0.029 0.023 0.013 0.008 0.004 0.003 

100 100 0.25 1.000 0.995 0.948 0.909 0.835 0.760 0.660 0.579 0.506 0.437 0.367 0.294 0.215 0.124 0.068 0.029 0.013 

10 100 0.025 1.000 0.989 0.886 0.807 0.643 0.494 0.367 0.302 0.251 0.206 0.167 0.128 0.089 0.045 0.019 0.005 0.002 

1 100 0.0025 1.000 0.982 0.860 0.758 0.535 0.282 0.124 0.099 0.081 0.065 0.051 0.040 0.028 0.013 0.006 0.002 0.001 

0.1 100 0.00025 0.999 0.985 0.885 0.827 0.707 0.552 0.293 0.069 0.012 0.009 0.008 0.007 0.006 0.004 0.003 0.003 0.002 



curve), with the corresponding bombardment solutions 
added as straight lines. Note that, contrary to what we 
found for T m ax, the first term in Eq. (|2^) does not repre- 
sent the bombardment solution, but rather the cyclotron- 
dominated shock heated plasma. It connects to the bom- 
bardment solution as the non-hydrodynamic regime is ap- 
proached and bridges a gap of two orders of magnitude in 
x s between the bombardment and bremsstrahlung solu- 
tions. Clearly, the quantitative determination of x s (m, B) 
requires radiation-hydrodynamical calculations. The in- 
dependence of Eq. ( |2(i| ) reproduces that of Eq. (|l(i|), the 
molecular weight dependence is added here and taken from 
Eq. (§. 



2.12. Geometrical shock height h s h 

Fig. ^ shows the quantity hgh B 2& for M w a = 0.6 M and 
for field strengths between 10 and 100 MG. We fit the data 

by 



h s \ 1 B 7 



2.6 



CQ/ic 



B. 



-2.6 



brcms 



1 1/7 



(21) 



where m as independent variable enters via ft s h, brcms 
from Eq. (|l4j) and the fit parameters Cq and 7 are 
listed in Table 2. The fits for 0.6 M Q and 1.0 M Q are 
shown (solid and dotted curve). The limiting dependen- 
cies for large and small rhB~ 26 , respectively, are h s b = 
h s h, brcms oc m _1 and /i s h= cq B^ 2 6 = const., as pre- 
dicted by Eq. (|ll|). The shock height is related to x s 
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Fig. 6. Column density x s of the post-shock cooling re- 
gion as a function of rh B^ 2e for the same parameters as 
in Fig. |^. The solid curve is the best fit from Eq. (poj). 
The dotted curve indicates the corresponding fit to the 
data for 1 M Q . The straight lines in the lower left denote 
the bombardment solutions from Eq. ( fi7| ) for 0.6 M and 
1.0 M Q , 

by the mean post-shock density p = x s /^sh which is 
~p~ = 6.97 po for the bremsstrahlung-dominated shock so- 
lution (see Eqs. [l3| and pT[ ). For the cyclotron-dominated 
shock-heated flow, the mean post-shock density increases 
to p = (bo/co)pov — 30 po, a result which can not be 
obtained from simple theory. 



3. Emission regions of finite lateral width 

In our radiation-hydrodynamical calculations, energy con- 
servation is enforced and the radiative luminosity of the 
infinite layer per unit area equals the accretion energy 
mv . Any real emision region, however, has a finite lateral 
width and looses energy not only from its top and bottom 
surfaces but also from its sides (Fig. 1). Two-dimensional 
radiation hydrodynamics would then be needed to calcu- 
late the temperature structure. In this section we discuss 
to what extent our one-dimensional results are still appli- 
cable to regions of finite extent. 



Table 2. Fit parameters of Eqs. (|19|), 
three values of the white dwarf mass M wc i 



and <ht) for 
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Fig. 7. Same as Fig. ^| but for geometrical shock height 
h a h contained in the quantity h^B 2 ® . The fits refer to 
0.6 Mq (solid curve) and 1 M Q (dashed curve). The straight 
lines represent the corresponding bremsstrahlung solu- 
tions of Eq. (ph. 
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Fig. 8. Overall spectral energy distributions for an emis- 
sion region on an 0.6M Q white dwarf at d = 10 pc 
with B = 30 MG, D = 10 8 cm, and m = 100,10,1, 
and 10 _1 gcm _2 s _1 (from top). The dashed sections in- 
dicate the emission directed away from the white dwarf 
at $ — 5° without the reflection albedo, the solid curves 
include the latter. The dotted curve for m = 1 gcm _2 s _1 , 
B = 100 MG indicates the increased cyclotron radiation 
and the reduced bremsstrahlung flux and temperature for 
this field strength. 

3.1. Emitted spectra 

We consider an emission region as depicted in Fig. 1, with 
finite width D, cross section D 2 , and field strength B. In a 
first step, we adopt the temperature and velocity profiles 
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along the flow lines, T(x) and v(x), calculated for infinite 
D also for the case of finite D. 

In the Rybicki code the radiation transport equation 
was solved with a mean cyclotron absorption coefficient 
and electron scattering was included. For the emission re- 
gion of finite extent, we calculate the outgoing flux at an- 
gle $ by ray tracing using the temperature profiles along 
slanted paths as shown schematically in Fig. 1. For rays 
starting or ending on the side surfaces, the temperature 
and density profiles were truncated appropriately. We ac- 
count separately for the cyclotron emissivities in the or- 
dinary ray (index o) and the extraordinary ray (index e) , 
and add 50% of the free-free emissivity with Gaunt fac- 
tor to both. We neglect electron scattering in the emis- 
sion region, but include the atmospheric albedo A v (van 
Teeseling et al. 1994). Each ray yields a contribution AI^ 
to the integrated intensity I„ (in ergs -1 Hz" 1 sr" 1 ) in that 
direction and the summation is extended over n rays, 

z„(0) - = E A(jW ( / e ^o^ s ) e- T - ( "' s) 

2=1 i—1 q 

+ J eW(#,s)e- T %(Wds), (22) 
o 

where r£*o(i?, s) and 7^*1 ($, s) are the optical depths of the 
ordinary and extraordinary rays along path i at angle i?, 
Smax is the pathlength along that ray and Aer^ the effec- 
tive projected area associated with it. The albedo contri- 
bution at i) is calculated as A v T v {v: — -d) and is not yet 
included in Eq. Ell. Reprocessing of the flux absorbed in 
the white dwarf atmosphere is not considered in this pa- 
per and the corresponding flux is, therefore, missing from 
our spectra. The spectral luminosity L v is obtained by in- 
tegrating Eq. ^ over An and the total luminosity L by 
integration over all frequencies. 

Fig. H shows the spectral flux at = 5° against the field 
direction emitted by an emission region with B = 30 MG 
and an area of 10 16 cm 2 (D — 10 s cm) on an 0.6 M@ white 
dwarf at a distance of 10 pc. Cyclotron emission domi- 
nates for low m and bremsstrahlung for high to. Free-free 
absorption becomes important near 10 15 Hz at the high- 
est to, but in reality this spectral region is dominated by 
the quasi-blackbody component produced by reprocess- 
ing of the incident flux in the white dwarf atmosphere. 
The results of WB96 on the ratio of the cyclotron vs. 
bremsstrahlung luminosities as a function of to and B 
remain basically valid, but will be modified if the shock 
is buried in the atmosphere and X-ray absorption is ac- 
counted for. 

Fig. ^| illustrates the optical depth dependence of the 
cyclotron spectra at d = 80°. Cyclotron emission lines at 
low to change into absorption features at high to. Since 
in real emission regions the fractional area of the high-m 
section is small (Rousseau et al. 1996) observed spectra 
show emission lines. 



io" 22 




v (Hz) 



Fig. 9. Cyclotron section of the spectral energy distribu- 
tions for the same set of parameters as in Fig. 8, except 
§ — 80°, and for mass flow densities to = 1,10 _1 , and 
10 _2 gcm _2 s _1 (from top). 

3.2. Specific luminosity 

An isolated emission region of lateral width D, shock 
height /i s h, and the temperature profile T(x) of the in- 
finite layer appropriate for the mass flow density to will 
have C = mv for optically thin and C > rhv for opti- 
cally thick emission. The overestimate in the latter case 
results from radiation emerging from the sides of the re- 
gion without a compensating influx. For the optical depths 
considered here, bremsstrahlung is practically free of such 
overestimate, cyclotron radiation is not. 

Let us assume for the infinite layer that mv feeds two 
components of C, namely Ahin and Ahick = mv — Ahin- 
For finite D, we then have 

£ ^ Ahin + (mv - Ahin) ^1 + '^j^J > mv , (23) 

where the second term in brackets is the ratio of the total 
surface area of the emission region over the sum of top 
and bottom areas. For simplicity, we have neglected the 
temperature variation over the surface of the emission re- 
gion and taken the energy loss per unit area as constant. 
We rearrange the terms in Eq. (|2^) to form a quantity A 
which relates C/mv to the aspect ratio of the emission 
region, h sil /D, 

A=(J--l)^L~l-^™. (24) 
\mv ) 2h sh mv a 

Fig. |l0| shows A as a function of mB^ 26 , calculated for 
model columns with h s h/D — 0.1,1 and 10 in the way 
described in the previous section. To a first approxima- 
tion, A is independent of h s h/D and the relative luminos- 
ity error (C — m)/rnv increases proportional to h^/D 
for a given to 73f 2 6 . The quantity A is negligibly small 
for large mB^ 2 ' 6 where bremsstrahlung dominates and 
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Fig. 10. Quantity A from Eq. ( |24| ) measuring the ex- 
cess luminosity C/mv as function of the aspect ratio of 
the emission region, hs^/D. The symbols refer to nar- 
row and pillbox-shaped emission regions with D/h s h= 0.1 
( ),D/h sh = 1( ),D/h sh = 10 ( ), the overplotted + and 
x refer to field strengths of 10 and 100 MG, respectively, 
the uncrossed symbols to B — 30 MG. 

reaches A ~ 0.4 for low mBj 2 ' 6 where cyclotron radia- 
tion dominates. Note that A never reaches the optically 
thick limit of unity because bremsstrahlung and optically 
thin cyclotron emission always contribute. To give an ex- 
ample, h^/D = 0.5 and A ~ 0.4 imply C/rhv ~ 1.4, i.e. 
an overestimate of the luminosity by 40%. 

In order to assess the size of the possible error for 
AM Her stars, we estimate h s h/D for a typical accretion 
rate of M = 10 16 gs _1 as a function of to. The lin- 
ear width of the emission region is D ~ (M/rh) 1 / 2 ~ 
lOVm 1 / 2 cm with to in gem s . For a brcmsstrahhmg- 
dominated flow, Eq. ( pT| ) yields /i s h,brcms/ D ~ 0.48/to 1 / 2 
(for M w( j = 0.6 Mq) which is less than unity since in this 
case to <; lgcm _2 s _1 . For a cyclotron-dominated flow, 
Eq. ( f2l| ) correspondingly yields /i s h,c yc /-D ~ 2m 1 / 2 /Bf 6 
which is again less than unity since now to ^ 1 gcm _2 s _1 
and typically By > 1. Hence, C is seriously overestimated 
only for isolated narrow subcolumns which are not radia- 
tively shielded by neighboring fluxtubes. 

3.3. Temperature correction 

Let the application of the unmodified one-dimensional 
temperature profile T(x) yield a specific luminosity C = 
r ■ mv with r > 1. We can then take then either: (i) iden- 
tify the parameters of this emission region with those ap- 
propriate for the increased mass flow rate ml — r ■ to; or 
(ii) recalculate the emission for a reduced mass flow rate 
to" = m/r and identify temperature and emission of that 
region as appropriate for the initial to. In case (ii), T max 
and x s are reduced to T max (m") and x s {rh"). This ap- 



proach demonstrates that the rising sections of the rela- 
tions displayed in Fig. | (Eq. |l|) and Fig. | (Eq. |o|) are 
further depressed for narrow columns, while the horizontal 
parts, where optically thin bremsstrahlung dominates, are 
not affected. Both approaches secure energy conservation 
but can not replace a proper treatment of the problem. 
They are not recommended for isolated tall columns. 

4. Discussion 

We have solved the equations of one-dimensional, two- 
fluid stationary radiation hydrodynamics for the shock- 
heated plasma in the emission regions on accreting mag- 
netic white dwarfs for a wide range of mass flow densities 
to = Pqv and field strengths B. For given B and M wc j, 
the peak electron temperature T max and the column den- 
sity x s of the emission region are physically related to 
to as the independent variable of the theory. They are 
no longer independent variables as in the frequently em- 
ployed "constant-A models". It is possible, therefore, to 
interpret the observed spectral energy distributions of ac- 
creting magnetic white dwarfs in terms of the distribution 
of mass flow densities present in their accretion spots. 

We now discuss to what extent the application of these 
results is limited by the simplifications made in our calcu- 
lations. One major simplification is the assumption of sta- 
tionarity which implies that we neglect the possible occur- 
rence of shock oscillations (Imamura et al. 1996, Saxton & 
Wu 1999, and references therein) and that we can not treat 
rapid time variability of rh. Since our approach can ac- 
commodate a range of m to occur in neighboring columns, 
the emitted spectrum will still approximate the true time- 
averaged spectrum if to varies only on time scales exceed- 
ing the post-shock cooling time r coo i ~ 4/i s h/u ~ Is. 
In the presence of shock oscillations which have periods 
of order r coo i, our results yield a mean temperature and 
column density which need not agree with the true time- 
averaged value if the oscillation is nonlinear (Imamura & 
Wolff 1990). 

The assumption of a one-dimensional flow implies that 
we neglect the convergence of the polar field lines of 
the white dwarf. In the spirit of this approximation, we 
have included the acceleration of the post-shock flow by 
a constant gravity g = GM w d/i? w d 2 and neglected the r- 
dependence of g. 

The assumption of an infinite layer implies that there 
is no temperature gradient perpendicular to the flow. This 
is no restriction for bremsstrahlung-dominated flows, but 
in such gradient is always established in columns of fi- 
nite width D by optically thick radiation components and 
lowers the mean electron temperature averaged across the 
column at any position x. We have suggested a simple 
first-order correction for the implied overestimate in T(x) 
which ensures conservation of energy and provides some 
remedy for narrow columns with /i s h ~ D. For very nar- 
row columns with h s h 3> D or absolutely tall columns with 
/i s h > 0.1 i?wd, the main radiative energy flow may be side- 
ways and the approach of Wu et al. (1994) becomes more 
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appropriate. In summary, our results are valid whenever 
h s h < -Rwd and h sh <, D. 

On the positive side, we consider our largely correct 
treatment of the two-fluid nature of the post-shock flow. 
One-fluid treatments (e.g. Chevalier & Imamura 1982, 
Wu et al. 1994) can account for cooling by cyclotron ra- 
diation in addition to bremsstrahlung, but are limited, 
by definition, to mass flow densities sufficiently high to 
ensure quick equilibration of electron and ion tempera- 
tures. They can not describe the substantial reduction of 
the peak electron temperature below the one-fluid value 
which we show to be present at low mass-flow density 
m and/or high magnetic field strength B. As a result, 
our description catches the essential properties of such 
columns: (i) dominant cyclotron cooling causes the peak 
electron temperature to stay far below the peak tempera- 
ture of the one-fluid approach; (ii) cyclotron cooling causes 
a drastic reduction in the column density and the geomet- 
rical shock height of the post-shock flow compared with 
pure bremsstrahlung cooling; and (hi) peak temperature 
and column density vary smoothly between the two limit- 
ing cases, the bremsstrahlung-dominatcd high-m regime 
(Aizu 1973) and the cyclotron-dominated low-m bom- 
bardment solution (WB92, WB93). The latter denotes the 
transition to the non-hydrodynamic regime and, gratify- 
ingly, our calculations recover the bombardment solution 
at the lowest mass flow densities accessible. Compared 
with WB96, we obtained numerically more accurate re- 
sults and have cast these into simple-to-use fit formulae 
which facilitate the modeling of emssion regions within 
the geometrical limitations noted above. No other two- 
fluid calculations with the full optically thick radiative 
transfer are available. 

The remaining, mainly geometrical limitations of our 
approach are inherently connected to the one-dimensional 
radiative transfer. Extension of the calculation to two 
dimensions encounters two problems: (i) a substantial 
increase in complexity and (ii) the introduction of an 
additional free parameter in form of the lateral width 
of the emission region. Therefore, we consider our one- 
dimensional approach with the correction explained in 
Sect. 3 as a reasonable compromise, with the noted ex- 
ception of tall columns. 

The present results can be used to quantitatively 
model the emission regions on accreting magnetic white 
dwarfs. The discussion of the overall spectral energy distri- 
bution of such objects requires to account for shocks being 
buried in the photosphere of the white dwarf. Such model 
allows to obtain the m— distribution in the accretion spot 
from observational data and is presented in Paper II of 
this series. The emission properties of AM Herculis bina- 
ries depend not only on m but vary also systematically 
with field strength: this dependence is described in Paper 
III. 
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